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Abstract 

The proton-inelastic scattering from ^'^He and ^'^^Li nuclei are studied in a folding model ap- 

o : 

^ . proach. A finite-range, momentum, density and isospin dependent nucleon-nucleon interaction 



(SBM) is folded with realistic density distributions of the above nuclei. The renormalization 
factors N/j and N/ on the real and volume imaginary part of the folded potentials are obtained 
^ . by analyzing the respective elastic scattering data and kept unaltered for the inelastic analysis 

at the same energy. The form factors are generated by taking derivatives of the folded poten- 
tials and therefore required renormalizations. The j3 values are extracted by fitting the p + 
^'^He/'^^Li inelastic angular distributions. The present analysis of p + '^He inelastic scattering 
to the 3.57 MeV excited state, including unpublished forward angle data (RIKEN) confirms L 
= 2 transition. Similar analysis of the p -|- ^He inelastic scattering angular distribution leading 
to the 1.8 MeV (L = 2) excited state fails to satisfactorily reproduce the data. 
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1. Introduction 

Recent advances in nuclear physics have given us the opportunity to delve into unique problems, 
hitherto unknown. An example is the neutron halo in the nucleus ^^Li, discovered as a consequence 
of its very large interaction radius, deduced from the measured interaction cross sections of ^^Li 
with various target nuclei 0, ^ ^. The halo of the nucleus extends its matter distribution to a 
large radius. Thus, the two valence neutrons in ^^Li, which form the halo, are extended well beyond 
the ^Li core and the 2n separation energy is exceedingly small (0.247 MeV). Besides the renowned 
example of ^^Li, there are other neutron rich nuclei, like ^He and '^He having extended valence neutron 
distributions called neutron halos/skins [0, ^. While ^He and ^'^Li have predominant core and two 
valence neutron structures, the ^He nucleus has four valence neutrons and has the largest neutron to 
proton ratio among these three nuclei. 

Along with the structures of such exotic nuclei near the drip lines, their excitation modes have also 
attracted considerable attention in recent times [§, ^ The existence of the neutron halo around 
the core, gave birth to the idea of a new resonance mode of excitation called, the soft dipole resonance 
(SDR) P, 1^, in which the halo neutrons oscillate against the core nucleus. From phenomenological 
and microscopic analyses of the ^^Li(p,p')^^Li^ 3 angular distribution data, excitation of the 1.3 MeV 
resonant state of ^^Li was found to correspond to L = 1 transition |T^, arising from SDR |Tl], 
Earlier microscopic analysis of p + ^Li inelastic scattering, include the work of Mani et al. |]13 



and Petrovich et al. where abnormally large deformations were predicted in [13|. For p + ®He 
inelastic scattering, coupled channel calculations using a Woods Saxon potential was carried out in 
Ref. [0 and spin-parity of the excited state (E* = 3.57 MeV) was predicted. The p + ^He inelastic 
angular distribution data leading to E* = 1.8 MeV ||T5| is also available. Since all of ^'^^Li and ^•'^He 
are loosely bound, their wave functions are quite extended in space. But, their difference in internal 
structures could lead to different excitation modes. In fact, SDR has been predicted in ^He and ^^Li, 
but not in ^Li and ^He. To appreciate this difference, a systematic study of their inelastic proton 
scattering data is desirable. 

In this work we present a consistent analysis of the proton inelastic scattering from ^'^^Li and ^'^He 
nuclei in a folding model approach. Although sophisticated folding calculations of two nucleon t- and 
g-matrices exist in a recent review of Amos et al. |jl6|, we feel that, a simpler calculation is always 
useful. The folding model is well known as a powerful tool for analyzing nucleus-nucleus scattering 
data at relatively low incident energies [jl4|, 0, |1^, |19[. It directly links the density profile of the 
nucleus with the scattering cross sections and is thus very appropriate for studying nuclei with exotic 
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matter distributions. However, in such analysis the choice of the nucleon-nucleon interaction is very 
crucial. As the unstable radioactive nuclei have different neutron and proton density distributions, 
an isospin sensitive nucleon-nucleon interaction is required to construct the folded potentials [1^, pO|] . 
For nuclei with low breakup thresholds, the folding model analysis also gives an estimation of the 
breakup channel coupling effects on the elastic through its requirement of a renormalization 
factor (N/j) on folded potentials. With increasing incident energy, the breakup channel coupling 
effect decreases, and the value of approaches 1.0. 

Using a density, momentum and isospin dependent finite range effective interaction |23] in a 
single folding model, the present work analyses the available low energy proton inelastic scattering 
data from ^Li (49.75A MeV) pl, ^^Li (68.4A MeV) [|], ^He (72.5A MeV) H, ||] and ^He (40.9A 
MeV) ||T^. A semi-microscopic analysis in the optical model (OM) framework is carried out for 
the p + ^Li,^He elastic scattering data while the OM analysis of p + ^^Li,^He elastic scattering 



at the above energies has already been performed on the same footing in the earlier work ||T9[. In 
the DWBA calculations of the nuclear excitation, with transferred angular momentum L, the form 
factors used are obtained by taking the derivative of the semi-microscopic potentials used. The p + 
^^Li inelastic scattering though analyzed in [1^, 11, 1^ is repeated here with microscopic real and 
volume imaginary potentials in addition to phenomenological surface and spin-orbit potentials, and 
by generating conventional form factors. Slightly different N/j and N/ values were obtained by a 
fit in [19| and these new values are used in the present work for the sake of completeness and 
comparison of the results with other nuclei. The formalism and the analysis are given in section 2 
while the summary and conclusions are given in section 3. 



2. Formalism and Analysis 

The form of the single folded potential |2^, used in the present work is, 
U(r1) = //9(r2)vNN(|ri - r2|)d^r2 

where, p{f2) is density of the nucleus and vnn is the effective interaction between two nucleons at 
the sites fl and with densities pi{ri) and P2('"l) respectively. A finite-range, density, momentum 
and isospin dependent effective interaction SBM (Modified Seyler Blanchard) is chosen, which has 
different strengths for pp (or nn) and pn interactions and its form is [p^ , 

Veff(r = |r-I -r2|,p,p) = -Q,,^[l - ^-d'ip. + p^n 

where, the subscripts '1' and 'u' refer to like-pair (nn or pp) and unlike-pair (np) interactions, respec- 
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tively. Here 'a' is the range of the two-body interaction in the configuration space, 'b' is a measure of 
the strength of repulsion with relative momentum 'p', while 'd' and 'n' are two parameters determin- 
ing the strength of density dependence. The parameters n, Ci, Cu, a, b, d are given in Table 1. These 



constants are found to reproduce the bulk properties of nuclear matter and of finite nuclei ^ 
and are known also to explain the p -|- 4,6,8jjg 6,7,9,iiLj scattering data successfully [l^, |lT], |12|, [T^, 
The parameters are determined without exchange effects and thus they contain the effect indirectly 
though in a very approximate way. 

The ^'^He and '''^^Li densities used in this work are shown in Fig. 1. The density prescriptions 
remain the same as that used in the earlier work on elastic proton scattering from these nuclei [|l^ 



For ^^Li, the cluster orbital shell model (COSM) density |jTO|, ^ has been used. The parametric 
form of the ^Li density is used from Ref. The density of ®He was also derived in the 

COSM approximation. It contains the extended distribution of valence nucleons and correspond to 
the experimental matter radius. For ®He, the p- inelastic scattering data at 40. 9A MeV is recently 
available [l^ and the L-transfer values are already predicted for some excited states |3^. In 
this work, a number of ground state densities derived [0, |32|, Q by using Faddeev wave function 
models called, PI, FC, FC6, Q3, Ql, FB, FA, K ,C, are employed to predict angular distributions for 
excitation to the 1.8 MeV state. Those density models incorporate different n-n and n-a potentials 
with a variation of the two-neutron separation energy E(2n) from about -1.15 MeV to -0.21 MeV 
and thereby a variation of the root mean squared (rms) radius of ^He. The rms radii corresponding 
to the above models are 2.32, 2.50, 2.53, 2.54, 2.56, 2.64, 2.64, 2.66, 2.76 fm respectively. These 
radii were computed assuming that the bare ^He core rms radius is 1.49 fm |3^. The ^He density 
is also plotted in Fig. 1 to show that all the nuclear densities have a tail extended well beyond the 
a-core. Since the interaction is isospin sensitive, separate neutron and proton densities of the nuclei 
are used |]T^, ^ ^ for folding calculations. 

Both the real (V) and volume imaginary (W) parts of the potentials (generated microscopically 
by folding model) are assumed to have the same shape, as in Ref. ||T9|, i.e. Vmicroij) = V + iW = 
(N/j -|- iN/)U(ri) where, N/j and N/ are the renormalization factors for the real and imaginary parts 
respectively |T^. These folded potentials with appropriate N/j and as required for elastic scat- 
tering fits (Table 3 of Ref. |]19| and this work), are used subsequently for inelastic scattering analysis 
in this work. The spin-orbit and the surface imaginary parts are taken from the phenomenological 
best fit calculations as before (Table 2 of Ref. [0 and this work). They needed minor adjustments 
in some cases for best fits as reported earlier |jl9| and in this work. The phenomenological potentials 
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had the following form, 

Vp^eno(r) = -Vo fo(r) - i f,(r) + 4 i a,W,(d/dr) f,(r) + 2{h/m^cy K.o 1/r (d/dr) f,.„(r) (L.S) + 

^ coul 1 

where, fx(r) = [1 + exp(^^^)]~^ and Rx = rxA^/^. The subscripts o,v,s,s.o denote real, volume 
imaginary, surface imaginary and spin-orbit respectively and Vq, W^, (W^) and V^.o are the strengths 
of the real, volume (surface) imaginary and spin-orbit potentials respectively. Vcoui is the Coulomb 
potential of a uniformly charged sphere of radius 1.40 K^/^. For each angular distribution studied be- 
fore ||19[ as well as here, best fits are obtained by minimizing x^/N, where = Y,^=i ""'''^^^ '{e^)^''^ ' 
<^th/ ciex are the theoretical/experimental cross sections at angle Ok-, AcXex is the experimental error 
and N is the number of data points. 

The best fit OM parameters for the p + ^Li elastic scattering at E = 49.75A MeV are given in 
Table 2. In the present semi-microscopic analysis, a search on N/j and N/ is carried out for minimum 
X^/N and the values are given in Table 3 and the fit is shown in Fig. 2a. The surface imaginary 
and spin-orbit potentials remain the same as obtained from the phenomenological best fits. For the 
p + ®He,^^Li elastic scattering data at 72. 5A and 68. 4A MeV respectively, the required N/j, N/ and 
phenomenological potentials (W^ and V^.o) are already available from the earlier analysis [1^. Using 
OM parameter setll of p + *He elastic scattering yields much lesser value and a better fit to 
the inelastic data. The corresponding best fit parameters for p + ®He scattering at 40. 9A MeV are 
also given in Table 2,3. To yield minimum x^/N, the r^ value had to be changed from 1.6 to 1.43, 1.32 
and 1.26 fm for the PI, Ql and C density models of ^He. These potentials are therefore used in the 
DWBA calculations of inelastic scattering with transferred angular momentum L. The calculations 



are performed by using the code DWUCK4 |3^. The conventional form factors, i.e. derivative of 
the potentials are used. The microscopic real and imaginary form factors have the same shape with 
strengths N^-'^ and Nf ^ respectively, where N^-^ = Nj:j^/rJ:^^, where the radius parameter is the 
rms radius of the folded potential. Thus the renormalization of the form factors is consistent with 
that for the folded potential. Form factors derived from phenomenological surface imaginary and 
spin-orbit potentials are also included. 

To fit the p + ^Li inelastic angular distributions leading to the 0.478 and 4.63 MeV excited states 
of ^Li, N/j_/ values from elastic scattering fits and the corresponding N^J"^ are employed. In the former, 
the best fit yields for the deformation parameter j3 a value of 0.59 (Fig. 2b, Table 3) for transferred 
angular momentum L = 2 (3/2^ to 1/2^). For the 4.63 MeV excited state of '^Li best fit gives a 
(5 = 0.82 for L = 2 (3/2^ to 7/2^). The calculations could explain the data up to 9cm ~ 110° and 



5 



value is calculated by incorporating the data points only up to that (Fig. 2c, Table 3). Earlier 
works 1^, |10|, 0, [1^ showed that the p + ^^Li inelastic scattering data at E = 68. 4A MeV could 
be satisfactorily explained for angular momentum transfer L = 1. The N/j = 0.50 and N/ = 0.15 
values are used in the present work as obtained from earlier work on p-elastic scattering |jl9[. The 
minimum test for best fit resulted in a (3 value of 0.58 (Fig. 3, Table 3), with L = 1. Changing the 
form factors generated from the surface imaginary and spin-orbit phenomenological potentials have 
negligible effects on the angular distribution. 

The ^He nucleus has a very high neutron-to-proton ratio. To confirm the spin-parity of the 
excited state E* = 3.57 MeV, studies over wide angular range are carried out. We try to fit both 
the reported ^He* data (setl) as well as the unpublished angular distribution data (setll) [p5|] . 
The setll data were obtained from invariant mass measurements as reported in |^. They were 
measured not by proton detection (like above), but by detecting ^He -|- n -|- n coincidences. Namely, 
the shape of the angular distribution for the inelastic scattering was extracted, while absolute value 
of the cross section was simply normalized in consistency with the above given inelastic data from 
proton measurements. It is seen that L = 2 (J'^ = 2+) gives best fit. L = 1 is excluded by the new 
measurement (setll) and L = 3 is excluded by all experimental points at larger angles (Fig. 4, Table 
3). The (3 = 0.28 is lesser than 0.44 in §. 

Similar analysis is extended to the recently available p -|- ^He inelastic angular distributions at 



40. 9A MeV (Fig. 5). The Faddeev wave function densities of He [plj pJj] are employed and the 
N/j, N/ values extracted from elastic scattering (Fig. 5a) are used. The angular distribution for the 
1.8 MeV (J'^ = 2"'") excitation is shown in Fig. 5b. For the sake of clarity only three calculations 
(PI, Ql, C models) are plotted and the rest lie between PI and C. It is found that contrary to very 
good fits to the inelastic data for p + ^He,'''^^Li, the present formalism is unable to give satisfactory 
fit to the p -|- ^He data. But the (3 value extracted from the optimized fit agrees closely with that 
of Aumann et al. and Rusek et al. |^8|] (/3 ~ 0.78). It is observed that though the three ^He 



densities correspond to a variation of rms radius from 2.32 fm to 2.76 fm, the corresponding change 
in the angular distributions of inelastic scattering is negligible. 

3. Summary and Conclusions 

A consistent folding model analysis of p-inelastic scattering of stable and unstable nuclei can 
provide valuable insight into their structure and reaction dynamics. The present work aims at 
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studying the inelastic scattering of protons from ^'^He and ^'^^Li nuclei. The semi-microscopic analysis 
followed here, involves a finite-range, momentum, density and isospin dependent nucleon-nucleon 
effective interaction (SBM) and realistic densities of different nuclei. Earlier, a similar analysis of 
p-elastic scattering data on these nuclei showed that renormalizations (N^j and N/) of the real and 
volume imaginary part of the folded potentials [|l^] are required. These factors, once determined from 
the elastic data, are kept unaltered for the inelastic data analysis to ensure a model independent 
study. 

The conventional way of generating the form factors is followed, i.e., by taking the derivatives of 
the potentials (microscopic real and imaginary as well as phenomenological surface imaginary and 
spin orbit). Deformation parameters are extracted from the analyses. The unpublished p5 



forward angle data of p + '^He inelastic angular distribution shows that best fit implies an L = 2 
transition (i.e J3 57 = 2+) whereas L = 1, 3 fail to reproduce the exact structure. 

Analysis on the same footing enables us to study elastic and inelastic p + ®He scattering at E = 
40. 9A MeV. In contrast to the other nuclei the p -|- ®He inelastic scattering to the 1.8 MeV state (J^g 
= 2+, L = 2) could not be satisfactorily explained by the present formalism, though the extracted /? 
value agrees closely with previous works |3^, Moreover, inclusion of various ground state density 
prescriptions of varying r.m.s radii (from 2.32 fm to 2.76 fm) have negligible influence on the inelastic 
observables studied here. 

The authors gratefully acknowledge A. A. Korsheninnikov for sending his unpublished experimen- 
tal data in a tabular form and giving his consent to use it in this paper. Thanks to I. J. Thompson 
for providing the ^He densities, F. Auger for sending the data in tabular form, R. Kanungo for the 
folding model code and S. K. Samaddar for useful discussions. D.G. acknowledges the CSIR, India 
for financial support. 
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Table 1: 

Parameters of the SBM interaction in MeV-fm units 



n 


Ci 




a b 


d 


2/3 


215.7 


669.3 


0.554 668.7 


0.813 
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Table 2: 

Optical potential parameters used in p + nucleus elastic scattering 



Nucleus 


E/A 




To 








(MeV) 


(MeV) 


(fm) 


(fm) 


(MeV) 


^Li 


49.8 


36.70 


1.210 


0.550 


5.62 


"Li 


68.4 


14.50 


1.385 


0.546 


4.26 


SRe 


72.5 


21.60 


0.970 


0.857 


1.43 


^He 


40.9 


45.40 


0.990 


0.612 


2.60 



(fm) 



a^, 

(fm) 



(MeV) 



(fm) 



(fm) 



V..O 

(MeV) 



(fm) 



(fm) 



J/A 
(MeV fm^) 



Ref. 



1.730 
0.560 
0.8 
1.101 



200 



1.220 
1.160 
633 
0.690 



4.60 
1.50 
3.47 



0.560 
1.600 
1.600 



1.160 
0.717 
0.772 



4.90 
5.90 
6.39 
5.90 



1. 

0.8' 
1. 
0. 



0000 
000 
4300 
6770 



.530 
.630 
.801 
.630 



424.6 
211.5 
246.1 
397.7 



m 
i 



\*] this work 



Figure Captions 

1. The densities of (a) ^'^•'^He (b) ^He, '''^^Li used in this work (see the text for references). 

2. The experimental angular distributions and folding model calculations (employing SBM inter- 
action) of ^Li at 49.75A MeV for (a) elastic and (b) E* = 0.478 MeV (1/2^), (c) E* = 4.63 
MeV (7/2~) state for inelastic scattering. Here L = 2 in (b) and (c). The corresponding N^, 
N/, N^"^, Nf^ values and phenomeno logical surface imaginary and spin-orbit parameters are 
given in Table 2, 3. 

3. The same as in Fig. 2 for ^^Li at 68. 4A MeV and E* = 1.3 MeV state for inelastic scattering. 
Here L = 1. 

4. The same as in Fig. 2 for ^He at 72.5A MeV for (a) elastic and (b) E* = 3.57 MeV state for 



inelastic scattering. In (b) the data represented by solid circles (Set 1) are from while the 
data represented by hollow circles (Set 11) are yet unpublished and obtained from pSf. The L 
= 1,2,3 calculations are shown by dashed, solid and dotted curves respectively. 

5. The same as in Fig. 2 for '^He at 40.9A MeV for (a) elastic and (b) E* = 1.8 MeV state for 
inelastic scattering. Here L = 2. The calculations are shown by the dashed, solid and dotted 
curves for the PI, Ql and C models. 
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Table 3: 

Renormalizations of SBM folded potentials and form factors for p-nucleus scattering at incident 
energy (E/A) and excited state energy (E*) in MeV, angular momentum transfer (L), deformation 
parameter (/?), volume integral (J/A) of the real folded potential in MeV fm^ and x^/N values from 
the elastic and inelastic scattering best-fits 



Nucleus 


E/A 


E* 


Ni^ 


N/ 


r^ 

rms 




Nf^ 


L 








J/A 


^Li* 


49.8 


0.478 


0.75 


0.24 


2.853 


2.140 


0.685 


2 


0.59 


5.681 


6.692 


423.8 


^Li* 


49.8 


4.630 


0.75 


0.24 


2.853 


2.140 


0.685 


2 


0.82 


5.681 


1.369 


423.8 


iiLi* 


68.4 


1.300 


0.50 


0.15 


3.909 


1.954 


0.586 


1 


0.58 


0.488 


0.353 


282.6 


8He* 


72.5 


3.570 


0.41 


0.00 


3.300 


1.353 


0.000 


1 


0.32 


0.594 


37.007 


247.7 


«He* 


72.5 


3.570 


0.41 


0.00 


3.300 


1.353 


0.000 


2 


0.28 


0.594 


0.452 


247.7 


«He* 


72.5 


3.570 


0.41 


0.00 


3.300 


1.353 


0.000 


3 


0.55 


0.594 


24.777 


247.7 


6He*(Pl) 


40.9 


1.800 


0.68 


0.00 


3.125 


2.125 


0.000 


2 


0.71 


6.247 


11.361 


444.1 


'He*(Ql) 


40.9 


1.800 


0.68 


0.00 


3.358 


2.283 


0.000 


2 


0.71 


6.907 


14.743 


469.8 


6He*(C) 


40.9 


1.800 


0.68 


0.00 


3.558 


2.419 


0.000 


2 


0.71 


7.980 


18.178 


484.8 
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